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Abstract 

A diagrammatic technique is derived for calculating processes involving the produc- 
tion of large numbers of particles. As an example, the amplitude of three particles 
scattering into n particles at the kinematical threshold in A</> 4 theory is calculated. 
From the 3 — > n amplitude, it is seen that the one-loop amplitude for 2 — > n processes 
is suppressed by a factor of 1/n. 



It has been shown that one can study multiscalar production at the tree-level by 
using recursion relationsE or the reduction formula methodfl . These techniques 
enable one to calculate the tree-level amplitude of one off-shell scalar particle decaying 
into many scalar particles at rest. 

These methods have been shown to be applicable to other problems. The one- 
loop correction was found in A0 4 theory for both the case of unbrokenS and broken@ 
reflection symmetry. In finding these corrections the amplitude of 2 on-shell scalar 
particles producing n particles at the kinematical threshold was discovered to vanish 
for n > 4 in the case of unbroken symmetry and n > 2 in the case of broken symmetry 
(for a discussion on this phenomena see ||). This effect was further discussed by 
Argyres, Kleiss, and PapadopoulosS and by Brown and ZhaiS. Argyres, Kleiss, 
and Papadopolous have discussed the cross channel of this process (i.e. one off- 
shell particle decaying into n particles at rest and one particle above the kinematic 
threshold). 

In this paper, we discuss a diagrammatic technique for calculating these ampli- 
tudes. We show how to build a set of Feynman rules for the generating functions 
of processes involving multiparticle production. With these rules, one can consider 
processes involving M scalars scattering into N scalars with none of the M or N 
particles neccesarily being at rest. However, instead of allowing one to calculate the 
M — > N amplitude as in normal diagrammatic techniques, this technique allows one 
to calculate the generating function for M — > N + n processes where the n particles 
are at rest. 

It must be emphasized that this procedure does not allow the study of many 
particles above the kinematical threshold. While large numbers of particles at rest 
may be produced, the calculations become more and more complicated as the number 
of particles above the threshold is increased. This approach does provide a compact 
method to calculate the amplitude for a few particles scattering to produce a few 
particles above the threshold and many particles at the kinematical threshold. 

We shall demonstrate this technique by calculating the 3 — > n amplitude in A0 4 
theory. Through the unitarity relations, this amplitude provides an easy way to study 
the one-loop properties of 2 — > n processes. 

We will be discussing A0 4 theory with the Lagrangian, 

C = \{d<pY-\m^-\<p\ (1) 
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We will primarily be deriving Feynman rules for the case of unbroken reflection sym- 
metry (m 2 > 0). Most of what is discussed may be trivially extended to the case of 
broken symmetry as well as theories dicussed in 0, [JTOj and ||. 

We wish to consider the amplitude for an interaction involving N particles with 
arbitrary four-momenta. During this interaction, n particles will be produced at 
the kinematic threshold. The particles produced at the threshold will have a four- 
momentum denoted as q = (1,0). The generating functions for these amplitudes can 
be denoted as the N-point Green function, G^\pi, ■■■,Pn)- To calculate this function, 
it is easiest to consider the amplitude of processes involving one particle decaying into 
N — 1 particles above the threshold and n particles at the threshold. This amplitude 
is written as 

4*V, ..., Pn ) = G^ +1 \-P - nq, Pl , ..., Pn ), (2) 

where P = p\ + ... + p n . 

The amplitude can be considered recursively as in Fig. 1. We can write the 
recursion relation 

A^n\pi, ■■■iPn) = (p \2~T E fin,ni+n 2 +n 3 E $N,Ni+N 2 +N 3 

[r + nq) i nun2in3 n 1} n 2 ,n 3 



„ .„ ,„ , E < 1 \{PnA)A^\{P n2 })Alf\{P n ^ (3) 
n-L<.n 2 in 3 \ perms . 

where {Pm} is a set of m of the original momenta and the second sum runs over all 
possible permutations of the momenta. 

Equation can be simplified by using the function, 

gf » (?1 ,.., p „),- i (|)" /2 ^W^) . (4 ) 

With this substitution, we can construct a differential equation for the generating 
function^, 



g (N \ Pl ,...,p N ;x) E^af)(p lr ,p^. (5) 

n=l 

It can be verified that for all N > 1, the generating function must satisfy, 



x 2 ^ + (2P + l)x-jj- + P 2 -1 



A E S N , Nl+N2+Ns g^\x)g^\x)g^\x), (6) 

Ni,N 2 ,N 3 
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where the momentum arguments have been omitted for simplicity and the sum is 
implied over all permutations of the momenta. 

The homogeneous part of equation (^j) contains g(°> . This coefficient is the gen- 
erating function for multiparticle production at the threshold, i.e. the amplitude for 
1 — > n processes as found in previous worksii. Using the fact that 

(n-l)/2 



<n|0(O)|O> = (2k + l)\ 



A 



(7) 



(/°) is found to be, 



9 



(0) 



X 



A (l-x\ 

To separate the homogeneous part of equation (H) , one must simply extract the terms 
in the sum with N\, N 2 , or N 3 = N. Equation @ can now be recast in the form of 
the homogeneous differential equation, 



d 



d 



x 2 — + (2P + l)x— + P 2 -1 



24x 2 



r 2\2 



gW(x) 



dx 2 v u ' ' dx (1 — .r- 

A E S NtNl+N2+N:i g^ (x)gW (x)g™ (x) . 

N U N 2 ,N 3 <N 



(9) 



As was noted by Voloshin in 0, the operator on the left hand side of equation 
can be restated in a familiar form by making the substitutions, 



x 
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{N) (p u ...,p N ;x) 



-Pot n ,N 



y (pi,-,PN;t). 



(10) 



This substitution transforms the homogeneous operator in equation @ into the 
Schrodinger operator for an exactly solvable potential. The generating function can 
now be written as the solution of 

6 



dt 2 



1 + 



y N {t) 



ii) 



cosh 2 (t)_ 

A £ 5 N>Nl+N2+N3 y N i(t)y N *(t)y N °(t), 

N!,N 2 ,N 3 <N 

where P 2 is square of the magnitude of the spatial part of P. 

The operator in equation ([12]) has been discussed in With cu = VP 2 + 1, the 
Green function for the operator in equation (|T^) was found to be 



G u (t,1?) 



2oo 



(12) 
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where 

= e^F(-2, 3; 1 ± W ; — (13) 

and -F(a, 6; c; a;) is the hypergeometric function. 

With the aid of this Green function, the solution of equation fll2l) is found to be, 

r+oo 

g W(t) = \e-* GUt,t') £ 5 N , Nl+N2+N3 y N ^t')y N2 (t')y N ^t')dt'. (14) 

N U N 2 ,N 3 <N 

Equation (|T4]) allows one to recursively calculate any g( N ' if one knows all g^ M ' 
for all M < N . Equation (|14]) can be summarized by a set of Feynman rules for 
calculating amplitudes with multiparticle production in the background. 

1. Draw all possible diagrams with the external legs labelled. Do not include any 
propagators that represent particles at rest. 

2. Assign a unique time label to every vertice and to the endpoint of every external 

leg. 

3. For every external leg, multiply the amplitude by a factor of e et , where e is the 
energy of the of the particle at the endpoint of the leg, and t is the time label at the 
endpoint of the leg. 

4. For every line, multiply by a factor of G u) {ti ) t2) where t\ and t<i are the time 
labels of the two endpoints of the line, and u is the expected on-shell energy of the 
momenta flowing through the line, i.e. u = \/k 2 + 1. 

5. For every vertice, multiply by a factor if A. In addition, for every three-point 
vertice multiply by a factor of g^°\t). This last factor is present to account for the 
production of particles at the vertice. 

6. Integrate over all time labels except for one. 

7. Multiply by a factor of e~^ eii , where t is the time label that was not integrated 
over and the sum includes the label at every external leg. 

8. Each diagram must be multiplied by an appropriate symmetry factor. 

Note that spatial momenta is conserved at every vertice and over all. Energy, 
however, is not conserved. This basic procedure can be used to construct Feynman 
rules in a number of theories including A0 4 in both the case of broken and unbroken 
symmetry. 

For the case of unbroken symmetry, the Green function for the propagator is given 
by equation (jl2[), and the vertice amplitude is given by equation (|]). The rules will 
generate g^ N '{t). To find the generating function of the amplitudes, one simply has 
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to make the substitution, 



z(t) ^/-e*. (15) 
A 



The amplitude is given by the n th derivative of g^ N '(z), 

Jn 

A^ip U ...,p N ) = —g( N \p u ...,p N ; Z ). (16) 

In theories with a broken symmetry, the propagator is given by®, 

qu t >) = KWtivnt - g) + fcwmwy - g ; (17) 

where 

/±(t) = e^F(-2, 3; 1 ± 2u; — (18) 

1 + e l 

The generating function for the amplitude of one particle production is given by, 

The amplitudes are found in the same manner as in the case of unbroken symmetry. 
However, in the case of broken symmetry, one must make the substitution, 

z(t) = -^|e*. (20) 

Rules for other theories may be constructed if one can calculate the propagator and 
the generating function for 1 — > n processes. 

As an example, we shall calculate the tree-level amplitude for threshold produc- 
tion of n particle by three on-shell particles in a theory with spontaneously broken 
symmetry. This amplitude is divergent when the incoming particles are on-mass-shell. 
To consider the amplitude of the real process, one must take the residue of the triple 
pole as the energy of all three incoming particles tends to its mass-shell value. The 
generating functional can be seen from the diagram in figure 2 to be, 

g V)(t) =6A J e £1 * i a i (t,t 1 )e e2t2 G W2 (t,t 2 )e^G Wl (t,t 3 )^^[^^ dt l( it 2 dt 3 . 

(21) 

Taking the on-shell limit corresponds to the limit — > Inspection of equation 
( |17|) reveals that in order to have poles in the amplitude, the time ordering of the three 
propagators is uniquely determined, t±, £2, £3 < t. Causality considerations would also 
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lead one to believe that this is the only time ordering capable of producing a physical 
amplitude. 

Equation (pT|) contains three integrals over hypergeometric functions. These are 
easiest to handle by writing the hypergeometric function in terms of the standard 
hypergeometric series. The integral of interest in equation Q2TJ) is 

jy^F(-2, 3 ; 1 - 2.; ^ = f_J~« ± ggj^*, (22) 

where (a) n is the Pochhammer symbol, 

. , T(a + n + 1) , . 

1 (a + 1) 

When the integrand in equation (22) is expanded in a series in e* , it is easy to 
see that only the first term in the series shall contribute to the pole as e — > —oo. The 
residue of the pole in equation (|22]) is the simple expression F(— 2, 3; 1 — 2a;; 1). We 
can now recast the integrals in equation fl2~T|) in the simple form, 

/oo I 
e e * G u (t, t')dt' = ———F(-2, -2 - 2a;; 1 - 2oo; -e 4 ). (24) 
-oo (1 + e t ) z 

Using the definition of z in equation (|20|), we can write an expression for the 
amplitude in equation (|2T|), 



9 <3»(,) = /A(l + ^) _ , /A A( 1+ ^) 

V2( 1 _ v ^/22)- V2 1 + 2^, 2ui ' 

_ J* + A < 1+ ^) (1 _ JIil±f5i, 4- A < 1 + ^) (25) 

V 2 l + 2w 2 2o> 2 A V2 1 + 2cj 3 2w 3 ; V ; 



Equation ( j2"5| ) is the generating function for the tree-level amplitudes of 3 — > n 
processes at the kinematical threshold. In the calculation of an amplitude, oo\, 0% 
003, and the number of particles are not four independent parameters. These four 
parameters can be reduced to two independent parameters through the application 
of momenta and energy conservation principles. 

The 3 — * n amplitudes contained in the generating function in equation (p5|) 
can be written in a compact form in the asymptotic limit, i.e. when 001,002,003, and 
n » 1. The tree-level amplitude for 3 — > n production at the kinematic threshold 
in the large n approximation in the theory with broken symmetry is 

2 (x\ {n+1)/2 n'n 6 

/(3 -"»^b (26) 
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Voloshin has noted the 2 — > n at the one-loop level can be related to the tree- 
level 3 — > n amplitudeQ. This reasoning is based on the argument that in the 
spontaneously broken theory, the imaginary part of 1 — ► n processes must vanish to 
all orders in perturbation theory. One can consider the unitary cuts on the one-loop 
1 — > n amplitude to derive the condition, 

21mA 1 (1 -> n) = J A 1 {I -> 2)A°(2 -> n)dr 2 + 
y A°(l -> 2)A 1 (2 -> n)rfr 2 + /" -> 3)A°(3 -> n)rfr 3 , (27) 

where the integration is done over the phase space of the intermediate particles, and 
the superscripts denote the number of loops in the amplitude. The first term on the 
right hand side of equation (|27|) is zero due to the vanishing of A (2 — > n) for all n > 2 
as noted in || . Since the left hand side of equation (|27| ) is zero, the one- loop 2 — > n 
amplitude may be related to the tree- level 3 — > n amplitude in equation (|25|). The 
integration over the phase space of equation (p5|) is non-trivial, but can determine 
whether the 2 — > n nullification noted in || will continue at the one-loop level. 

The expression in equation (|26|) can now be combined with the 1 — >• 3 amplitude 
to find the contribution of the second term in the unitary relation, equation (p7|). This 
contribution is of the order of \^ n+l ^ 2 n\n 2 . The order of magnitude of the one-loop 
correction to the 2 — > n amplitude is, therefore, 

A\2-±n) = A (n+2)/2 n!n 2 . (28) 

Although the tree-level contribution vanishes in A0 4 theory, the one-loop contribution 
can be compared with the tree-level amplitude for 2\ — > n<p in an interactive theory 
as calculated by Brown and Zhai@. They found that the tree level 2 — > n amplitude 
should grow like X n ^ 2 n\n. 

The above arguments indicate that although the 2 — > n tree-level threshold am- 
plitude may vanish, diagrams at the one-loop level provide a correction of order Xn. 
If the loop expansion were a series expansion in powers of Xn 2 , the results of || would 
indicate that the one-loop 2 — > n correction should grow like A^ n+1 ^ 2 n!n 3 . There is a 
1/n suppression in the one-loop correction to threshold 2 — > n processes. 

The rules discussed in this paper provide a compact method for calculating many 
processes involving large numbers of particles at the kinematic threshold. They pro- 
vide a useful tool for investigating suppression of amplitudes in multiparticle pro- 
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cesses. They may prove to be of further assistance in studying multiparticle phenom- 
ena such as the vanishing of 2 — ■> n amplitudes. 

I would like to thank M.B. Voloshin for several helpful discussions. 
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Figure 1 - Symbolic representation of recursion relation 
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Figure 2 - Feynman diagram for 3 — > n processes 
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